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The statistics of domain walls for ground states of the 2D Ising spin glass with +1 and -1 bonds 
are studied for L x L square lattices with L < 20, and x = 0.25 and 0.5, where x is the fraction of 
negative bonds, using periodic and/or antiperiodic boundary conditions. Under these conditions, 
almost all domain walls have an energy Edw equal to or 4. The probability distribution of the 
entropy, Sdw, is found to depend strongly on Edw The results for Sdw when Edw = 4 agree with 
the prediction of the droplet model. Our results for Sdw when Edw = agree with those of Saul 
and Kardar. In addition, we find that the distributions do not appear to be Gaussian in that case. 
The special role of Edw = domain walls is discussed, and the discrepancy between the prediction 
of Amoruso, Hartmann, Hastings and Moore and the result of Saul and Kardar is explained. 



I. INTRODUCTION 

There continues to be a controversy about the na- 
ture of the Ising spin glass. The Sherrington-Kirkpatrick 
modelji with its infinite-range interactions between the 
spins, is described by the Parisi replica-symmetry break- 
ing mean-field theory^ 3 To understand models with 
short-range interactions on finite-dimensional lattices, 
however, it is necessary to include the effects of interfaces, 
which do not exist in a well-defined way in an infinite- 
range model. The droplet model of Fisher and Husep*^ 
which starts from the domain- wall renormalization group 
ideas of McMillanLSiS and Bray and Mooreji2*ii and 
studies the properties of interfaces, provides a very dif- 
ferent viewpoint on the spin-glass phase. More recently, 
numerical resultsi2ii&±i indicate that the actual situation 
for three-dimensional (3D) models combines elements of 
both pictures in a nontrivial way. 

In two dimensions (2D), the spin-glass phase is not sta- 
ble at finite temperature. Because of this, it is necessary 
to treat cases with continuous distributions of energies 
(CDE) and cases with quantized distributions of ener- 
gies (QDE) separately. 11 i 15 

In three or more space dimensions, where a spin-glass 
phase is believed to occur at finite temperature, the 
general framework of thermodynamics requires that the 
CDE and the QDE should be treated on the same foot- 
ing. The way this comes about is that in these cases 
the typical domain wall energy increases as a positive 
power of the size of the lattice. Thus the quantization 
energy becomes a negligible fraction of the domain wall 
energy for large lattices. All bond distributions behave 
in a qualitatively similar way, except that the QDE have 
finite ground state entropies &ii 

Amoruso, Hartmann, Hastings and Moored have re- 
cently proposed that in 2D there is a relation 

^ = 1+ 4(4)' (1) 

where df is the fractal dimension of domain walls, and 
0e is the exponent which characterizes the scaling of the 



domain wall energy with size. For continuous energy dis- 
tributions the existing numerical estimates of df and 9e 
satisfy Eqn. (1). 

It is known from the droplet theory^ that for the QDE, 
which have a positive entropy at zero temperature, 

d f = 29 s . (2) 

6s is the exponent for the scaling of domain wall entropy 
with size. Thus, for the QDE, Eqn. (1) provides a rela- 
tion between the scaling of the energy and the entropy 
of domain walls. It is not known how to calculate df 
directly for the QDE case, so we need to use Eqn. (2) to 
verify Eqn. (1) in that case. 

For the QDE, it is known that 6 E = OiiSiii Then using 
Eqn. (1) gives df = 5/4, or using Eqn. (2), S = 5/8. 
The calculation of 6s by Saul and Kardarji2*i2i found 
6s = 0.49 ± 0.02. Since df cannot be less than 1, this re- 
sult was interpreted as a strong indication that 6s = 1/2. 

In this work we will demonstrate that Eqn. (1) actu- 
ally does work for both continuous and quantized distri- 
butions in 2D. The actual behavior of the QDE probabil- 
ity distributions under finite-size scaling turns out to be 
more subtle than what has been assumed until recently. 20 
As pointed out by Wang, Harrington and Preskill^Si do- 
main walls of zero energy which cross the entire sample 
play a special role when the energy is quantized. 

We will analyze data for the domain wall energy, Ed w , 
and the domain wall entropy, Sdw, for the ground states 
(GS) of 2D Ising spin glasses obtained using methods 
from earlier workiiSii^iSS* 2 ^ We will demonstrate that 
for L x L square lattices the Edwards- Anderson 2 ^ (EA) 
model with a ± J bond distribution has a strong corre- 
lation between Edw and Sdw for the GS domain walls. 
Because of this correlation, we will need to treat domain 
walls of different energies as distinct classes, whose en- 
tropies scale in different ways as L is increased. We will 
find that the scaling parameter identified by Saul and 
Karda»i2ii£ is the one associated with domain walls hav- 
ing Edw = 0. It is not, however, the one which controls 
the dominant behavior for large L. 
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II. THE MODEL 

The Hamiltonian of the EA model for Ising spins is 

H = - JijSiSj (3) 

(ij) 

where each spin S% is a dynamical variable which has 
two allowed states, +1 and -1. The (ij) indicates a sum 
over nearest neighbors on a simple square lattice of size 
L x L. We choose each bond to be an independent 
identically distributed quenched random variable, with 
the probability distribution 

P(J l0 ) = x8(J lJ + 1) + (1 - x)8(J lJ - 1) , (4) 

so that we actually set J = 1, as usual. Thus x is the 
concentration of antiferromagnetic bonds, and (1 — x) is 
the concentration of ferromagnetic bonds. 

The data analyzed here used an ensemble in which, for 
a given value of x, every Lx L random lattice sample had 
exactly (1 — x)L 2 positive bonds and xL 2 negative bonds. 
Details of the methods used to calculate the GS energies 
and the numbers of GS have been described earlier^ and 
the data analyzed here were used before to obtain other 
properties of the models 



III. GROUND STATE DOMAIN WALLS 

The GS entropy is defined as the natural logarithm of 
the number of ground states. For each sample the GS 
energy and GS entropy were calculated for the four com- 
binations of periodic (P) and antiperiodic (A) toroidal 
boundary conditions along each of the two axes of the 
square latticed We will refer to these as PP, PA, AP 
and AA. In the spin-glass region of the phase diagram, 
the variation of the sample properties for changes of the 
boundary conditions is small compared to the variation 
between different samples of the same size^ except when 
x is close to the ferromagnetic phase boundary and the 
ferromagnetic correlation length becomes comparable to 
L. 

We define domain walls for the spin glass as it was done 
in the seminal work of McMillan. 8 We look at differences 
between two samples with the same set of bonds, and the 
same boundary conditions in one direction, but different 
boundary conditions in the other direction. Thus, for 
each set of bonds we obtain domain wall data from the 
four pairs (PP,PA), (PP,AP), (AA,PA) and (AA,AP). 
The reader should remember that the term "domain 
wall", as used in this work, refers only to this procedure. 
Saul and Kardaipi&IS follow the same procedure used in 
this work, but use the term "defect" instead of "domain 
wall" . 

It is important to realize that the meaning of a domain 
wall is very different when the GS entropy is positive, as 
in the model we study here, as compared to the standard 



TABLE I: Domain wall energy statistics for x = 0.25. Naf 
is the number of negative bonds in each sample, and N sam is 
the number of distinct bond configurations studied for a given 
L and Naf- Tit is the number of domain walls of each type 
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case of a doubly degenerate ground state. In the stan- 
dard case one can identify a line of bonds which forms a 
boundary between regions of spins belonging to the two 
different ground states. It is not possible to do that when 
there are many ground states. Despite this, we continue 
to use the term "domain wall" . 

Due to the P(Jij) we have chosen and the fact that we 
have only used even values of L, the energy is always a 
multiple of 4. Thus, the energy difference, Ed w , for any 
pair must also be a multiple of 4. The sign of Ed w for 
a pair is essentially arbitrary for the values of x we will 
study here, which are deep in the spin-glass region of the 
phase diagram. Thus we can, without loss of generality, 
choose all of the domain-wall energies to be non-negative. 
For x = 0.25 and x = 0.5, it turns out, crudely speaking, 
that about 75% of the time we find Edw — 0, and 25% of 
the time Edw = 4. For a given L, the fraction of Edw = 
domain walls, fo, is about 3% higher at x = 0.5 than at 
x = 0.25. It is interesting to note that Wang, Harrington 
and PreskillSi use an analytical argument to predict that 
fo is approximately 0.75, independent of x, in the spin- 
glass regime. The value of fo tends to increase slightly 
as L increases, but our statistical errors are too large for 
a good estimate of the finite-size scaling behavior to be 
made. Our results for x = 0.5 are consistent with the 
results of Amoruso et alm^ 

Our calculated statistics for these values of x 7 as a 
function of L, are given in Table I and Table II. Estimat- 
ing the statistical uncertainties in these numbers is not 
trivial, due to the fact that the values of Ed w obtained 
from the same set of bonds with the four different pairs 
of boundary conditions are not statistically independent. 
An upper bound on the statistical uncertainties is ob- 
tained by counting the number of samples, rather than 
the number of pairs. 

No domain walls with energies greater than 8 were ob- 
served at any L for these values of x. This, however, does 
not have much fundamental significance. The probability 
distribution for Ed w is also a function of the aspect ratio 
of the latticed 

The domain wall entropy, Sd w , is defined, by analogy 
to Edw, to be the difference in the GS entropy when the 
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TABLE II: Domain wall energy statistics for x — 0.5. Column 
labels as in Table I. 
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boundary condition is changed along one direction from 
P to A, with the boundary condition in the other direc- 
tion remaining fixed. The probability distribution of Sdw 
for the cases where Ed w = should be symmetric about 
0, and our statistics are consistent with this. As shown by 
Saul and Kardarji&iSi the variance of Sdw when Ed w = 
increases with L in approximately a linear fashion. This 
is illustrated in Fig. 1, for both x = 0.25 and x = 0.5. 
The calculated skewness of these essentially symmetric 
distributions is, naturally, consistent with zero, but their 
kurtosis is not. For x = 0.25 the calculated kurtosis for 
different values of L varies between 1.0 and 2.5, with no 
apparent trend in the L dependence. For x = 0.5, where 
the number of samples is about twice as large, the kur- 
tosis varies between 1.35 and 2.05, again without any 
apparent trend with L. If we average over L, we obtain 
1.58±0.18 for x = 0.25 and 1.70±0.10 for x = 0.5. From 
this we conclude that the Ed w = distributions are not 
Gaussian, and that they probably do not become Gaus- 
sian even as L becomes large. The results are consistent 
with a probability distribution which has a kurtosis of 
about 1.7, independent of L and x. Since the model is 
critical at T = for both values of x, this result is not 
surprising. 

When Edw is not zero, the relative signs of Ed w and 
Sdw are not arbitrary. Having chosen Edw to be nonncg- 
ative, we then find that, when Edw is positive, it turns 
out that Sdw is almost always positive. In Fig. 2 we 
show the behavior of the average value of Sdw for the 
cases where Ed w = 4, as a function of L. We see that 
for Edw = 4, the average value of Sd w (L) grows approx- 
imately as L 63 , for both values of x. More precisely, a 
least-squares fit to the form 



gives 



for x — 0.25, and 



Sdw(L) = AL e 



9 S = 0.639 ± 0.036 



0.628 ±0.027 



(5) 



(6) 



(7) 



for x — 0.5. These numbers should be compared to the 
prediction of droplet theory, 6 Eqn. (2). Because of the 
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FIG. 1: Variance of the entropy vs. L for the Edw ~ 
domain walls, for the cases x = 0.25 and x — 0.5. 



E=4 Domain Walls 
Avg Entropy vs. L 
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FIG. 2: Average entropy vs. L for the Edw = 4 domain 
walls, for the cases x = 0.25 and x = 0.5, log-log plot. The 
error bars indicate 1 a, and the statistical errors are larger 
for x = 0.25 because the numbers of samples are smaller. 



large GS degeneracy in the ± J Ising spin glass, one does 
not know how to compute df directly for this model. 

For the Ising spin glass with a Gaussian bond distribu- 
tion, however, this calculation in 2D is straightforward, 
and the result 26 i 27 i 28 is df = 1.27 ± 0.01. It would be 
desirable to improve the accuracy of our estimate of 9g 
sufficiently to demonstrate that 9s is really different for 
the CDE and the QDE, as predicted by Eqn. (1). The 
predicted numerical difference is small, however, so this 
will be difficult. 

The skewness and kurtosis of the Sdw distributions for 
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E=4 Domain Walls 
<S>/sd(S) vs. L 
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FIG. 3: Average entropy divided by its standard deviation 
vs. L for the Edw = 4 domain walls, for the cases x — 0.25 
and x = 0.5. 



Edw — 4 are both small, in contrast to the Edw = case. 
It is likely that these distributions become Gaussian in 
the large L limit. 

The result of Saul and Kardar ji&iii obtained by looking 
at the distribution of Sdw for all values of Ed w combined, 
was 9s — 0.49 ± 0.02. To obtain this exponent, Saul and 
Kardar fit the data at small values of Sdw For practical 
purposes, this part of the data belongs to the Ed w = 
component. It is clear that if our extrapolations from 
small L are correct, however, eventually, at large enough 
L, the width of the total Sdw distribution will be dom- 
inated by the Ed w = 4 component. Although they did 
not identify it, the peak in the probability distribution of 
Sdw due to the Ed w = 4 component is clearly visible in 
Fig. 19 of the 1994 paper of Saul and Kardar ~- 

In Fig. 3 we show the values of the ratio of average 
value of the S(L) distribution for Ed w = 4 divided by 
the rms width of the distribution. Except for L = 6, 
where the ratios are smaller, we see that this ratio is 
about 2.9 for x = 0.25, and about 3.2 for x — 0.5. Thus 
this characteristic of the distributions seems to depend 
on x. However, it seems likely that the ratio is actu- 
ally L-dependent, and that the large L limit is universal, 
independent of x. 



IV. DISCUSSION 

What we have learned is that in this model there are 
two distinct classes of domain walls, the Edw — domain 
walls and the Edw > domain walls. For x = 0.25 or 
x = 0.5, the only type of Edw > domain walls which 
we need to worry about are the Ed w = 4 ones. If we 
wanted to study values of x of 0.15 or less,- 9 or values 
of the aspect ratio less than one^£ then we would need 



to consider larger values of Edw It should be expected 
that these are qualitatively similar to the Ed w — 4 case, 
as explained below. 

As we have seen, the Ed w — 4 domain walls behave in 
a way which is essentially consistent with the predictions 
of the droplet model, but the Ed w = domain walls do 
not. This difference in behavior cannot be understood 
within the droplet model. It must be related to the spe- 
cial feature of the model caused by the quantization of 
the energy spectrum, which is the positive GS entropy. 
Another fact which we need to use is the recently discov- 
ered result that the average GS entropy is an increasing 
function of the GS energy. 

For an Edw > domain wall, a contribution to Sdw 
comes from this shift in the average GS entropy with the 
shift in the GS energy. What remains to be understood 
is why Sdw should scale with increasing L in the way pre- 
dicted by the droplet model. The conventional derivation 
of the droplet model^ appears to use the assumption that 
the GS is unique, up to a reversal of the entire state, in 
an essential way. The extension of the droplet model to 
the more general case was given by Fisher and Huse£ 

As long as Edw > 0, the two boundary conditions 
which we are comparing are not on an equal footing. As 
Wang, Harrington and Preskill 2 ^ express the situation, 
the Edw > domain wall does not destroy the topologi- 
cal long-range order. However, in the Ed w = case the 
two boundary conditions are on an equal footing, and the 
topological order is destroyed. Therefore the Ed w = 
class of domain walls can be expected to behave in a spe- 
cial way, which differs from the prediction of the droplet 
model. 

It is well established that this model displays large cor- 
rections to finite-size scaling^ The reader may object 
that nothing reliable about the large L behavior can be 
learned from data which is restricted to L < 20. How- 
ever, our results indicate that the reason for the large 
corrections to scaling is that the differences in the scaling 
exponents for different classes of domain walls is small. 
When each class is analyzed separately, this effect is elim- 
inated. Thus, going out to only L — 20 is not as bad as 
it might appear at first. A more significant limitation in 
our work is that the number of samples we have data for 
is much smaller than we would like. 

We also wish to point out that there is a natural simi- 
larity between the Ed w — domain walls in the 2D Ising 
spin glass with ±J bonds and the large-scale low en- 
ergy excitations which have been founcU2*A£ in the three- 
dimensional version of the same model. A similar, and 
possibly related large-scale low energy excitation has also 
been seen in the 3D random-anisotropy XY model. 32 
Similar ideas have also been discussed by Hatano and 
Gubernatisii^ 

The domain-wall renormalization group of McMillan 7 
is based on the idea that we are studying an effective cou- 
pling constant which is changing with L. For the case of 
a continuous distribution of bonds, we can use the energy 
as the coupling constant. For the quantized energy case, 
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what we should really do is a slight generalization of this 
idea. We should think of the coupling constant as the 
free energy at some infinitesimal temperature. When we 
do this, the entropy contributes to the coupling constant. 
And as we have seen, the distribution of Edw rapidly be- 
comes essentially independent of L as L becomes large. 
Under these conditions, it becomes natural to treat each 
value of Ed w as a separate class, representing a different 
coupling constant. 

Therefore, since Sdw increases as a positive power of L 
for Ed w — 4, and presumably for all other (finite) pos- 
itive values, this coupling constant must eventually be 
controlled by Sdw And this coupling constant is also the 
one described by the droplet model. 

The Edw = class behaves differently, and therefore 
represents a different coupling constant. It is natural to 
wonder if topological long-range order can be related to 
replica-symmetry breaking, and if the Ed w = domain 
walls can be described by the replica-symmetry breaking 
theory. We will not attempt to do this here. 

The argument that all classes other than the Edw = 
class behave in the same way should certainly be checked. 
The way to do this would be to do calculations with odd 
values of Li^ In this way we would obtain results for the 
Edw = 2 class, and possibly also for the Ed w = 6 class. 

The computing effort used to obtain the results de- 
scribed here was rather modest by current standards. 
The algorithm^ is capable of generating data for greater 
numbers of samples out to somewhat larger values of 
L (perhaps to L = 32) with a reasonable effort. The 
method of Lukic et alm& appears to be more powerful, 



and has been used by them to generate energy-entropy 
statistics for ground states of this model out to L = 50. 



V. SUMMARY 

We have studied the statistics of domain walls for 
ground states of the 2D Ising spin glass with +1 and 
-1 bonds for L x L square lattices with L < 20, and x = 
0.25 and 0.5, where x is the fraction of negative bonds, 
using periodic and/or antiperiodic boundary conditions. 
Under these conditions, almost all domain walls have an 
energy Ed w equal to or 4. The probability distribution 
of the entropy Sdw is found to depend strongly on Edw 
The results for Sdw when Ed w = 4 agree with the predic- 
tion of Amoruso, Hartmann, Hastings and Moore, 16 Eqn. 
(1). Our results for Sdw when Ed w = agree with those 
of Saul and Kardar^iS but in addition we find that the 
distributions are not Gaussian in that case, even in the 
limit of large L. Due to the special role of the Edw = 
domain walls, we can understand the difference between 
the scaling exponent found by Saul and Kardar and the 
prediction of Eqn. (1). 
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